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ABSTRACT 

We present the updated orbital solution for the transient Be X-ray binary V 0332-1-53 complementing historical measurements with the 
data from the gamma-ray burst monitor onboard Fermi obtained during the outburst in June-October 2015. We model the observed 
changes in the spin-frequency of the pulsar and deduce the orbital parameters of the system. We significantly improve existing constrains 
and show that contrary to the previous findings no change in orbital parameters is required to explain the spin evolution of the source 
during the outbursts in 1983, 2005 and 2015. The reconstructed intrinsic spin-up of the neutron star during the latest outburst is found 
to be comparable with previosly observed values and predictions of the accretion torque theory. 

Key words, pulsars: individual: - stars: neutron - stars: binaries 


1. Introduction 

The bright X-ray transient source V 0332-153 was discovered 


by the Vela 5B satellite during an outburst in 1973 (Terrell & 
|Priedhorsky|1984 1 . Observations with EXOSAT ten years later 


allowed the position of the source to be determined the X-ray 
pulsations with period of ~ 4.4 s ( Stella et al. [1985) 1 modulated 
by motion in an eccentric orbit (e ~ 0.3) with period of ~ 34.3 d 
to be detected. The optical counterpart of the X-ray pulsar has 
been identified to be a Be star (Honeycutt & Schlegel|1985 1 at a 


distance of ~ 7 kpc (|Negueruela et al.|1999[l. 


Typically for Be systems, V 0332-153 exhibits two types of 
ray outbursts with peak luminosities of ~ 10^^ erg s ' and 
10^^ erg s ' associated with enhanced accretion onto the neu- 


ergs 

tron star from the Be circumstellar disk close to the periastron. 
The so-called “normal” or Type I outburst occur regularly as the 
neutron star passes throngh the Be disk and have comparatively 
low luminosity. Longer and more luminous type II or “giant” out¬ 
bursts are more rare and thought to be related to the formation of 
a stable accretion disk around the neutron star during a periastron 
passage (Okazaki et al.|2002{ Martin et al.|2014). Eonr type II 


outbursts and a number of type I outbursts have been detected 


from the source so far (Terrell & Priedhorsky 

1984 Tsunemietal. 

1989 

Swank et al.|20041|Camero-Arranz et a 

.120151 |Doro shenko | 

et al. 

2015 Nakajima et al.|2015|l. 


In this note we report on the analysis of timing properties of 
V 0332-153 during the most recent type II ontburst which took 
place in June-October 2015 together with older observations and 
provide the updated orbital solution for the system. 


2. Data analysis and results 


Orbital parameters of V 0332-153 were first determined by [Stella 
jet al.|(|1985|) using th e EXOSAT data. Later jZhang et al.| ( |2005| l 
and jRaichur & Paul ( |2010| l used Rossi X-Ray Timing Explorer 
(RXTE) and International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) observations performed during the 2004-2005 gi¬ 


ant outburst to refine the orbital solution. These authors found 
significantly larger projected semimajor axis a sin i value than 
reported by Stella et al. ( [1985 I which was associated with the 
apsidal motion in the system. 


The most recent outburst in 2015 has been monitored with 
several instruments. However, the best timing information is 
provided by the Gamma-ray burst monitor (GBM) on board Fermi 
( jMeegan et al.|2009| l which regularly measured spin frequency 
of the source every 1-3 days throghout the outburst. The GBM 
consists of 12 Nal and two BGO detectors providing the full-sky 
coverage in 8 keV-40MeV energy range and is designed to 
detect and localize gamma ray bursts. However, it also proved 
to be very useful for detection and monitoring of pulsed signals 
from X-r^ pulsars. In particular, the Fermi GBM team regularly 
publishe^spin histories of selected pulsars including V 0332-1-53 
based on the analysis of CTIME data from Nal detectors in 8- 
50keV energy range. 


Pnlsations from V 0332-1-53 were first detected by GBM 
on MJD 57194 and the data used in this publication covers the 
interval until MJD 57290, i.e. almost three full orbital cycles. 
The orbital modulation of the spin frequency superimposed by 
a spin-up trend is clearly visible in the raw data shown in Eig. 
We use these measurements to constrain the orbital parameters 
and intrinsic spin evolution of the pulsar using the same approach 
as [Zhang et al. (2005 i. In addition, we repeat the analysis of 
RXTE data carried out by Raichur & Paul ( [20I0[ l. In particular, 
we reconstruct the spin history of the pulsar during the giant 
outburst in 2005 using the epoch-folding period search and RXTE 
PGA lightcurves in 3-21 keV energy range between MJD 53332 
and 53432. Using the obtained spin frequency measurements 
we followed procedures described in Zhang et al. (2005 I and 
[Raichur & Paul P010[ l to constrain the orbital parameters of the 
system. We modeled simultaneously the data from 2005 and 2015 
outbursts together with historical pulse frequency measurements 


' http://gammaray.msfc.nasa.gov/gbm/science/pulsars.html 
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Fig. 1. Observed pulse frequency modulated by orbital motion as measured with (left to right) Fermi GBM, RXTE and EXOSAT and 
Tenma during the three major outbursts from the source. Reconstructed intrinsic pulsar frequency (black dashed line) modulated by 
motion along the orbit with best-fit parameters (red line) together with fit residuals are also shown. 


reported by|Stella et al.|(|1985|l and|Makishima et al.|(|1990|) for 
the outburst in 1983-1984. 

We emphasize that both the orbital modulation and intrinsic 
variation of the pulsars spin frequency must be modeled in order 
to constrain the orbital parameters of the system. The intrinsic 
spin evolution of V 0332+53 during an outburst is rather compli¬ 
cated as found by |Zhang et al.| ( [2005) ; |Raichur & PauT| ( |2010| l who 
had to include the second pulse period derivative to describe it 
adequately. This is to be expected as the accretion torque spinning 
up the neutron star depends on the accretion rate which changes 
dramatically and also true for the latest outburst. In fact, we were 
unable to describe the spin history measured with Fermi GBM 
even with inclusion of the third spin frequency derivative as a 
parameter. 

Therefore, we opted to model the intrinsic spin evo¬ 
lution of the pulsar as a smooth interpolating func¬ 
tion (i.e. using the piecewise cubic hermite interpolating 
polynomial) defined by frequency values at fixed times 
Ti = M JD57194,57208,572 20,57230,57250,57270,57290. 
Following |Zhang et al.|(|20051l we included a systematic uncer¬ 
tainty of 7 X 10 “ Hz and 7 x 10 ^’ Hz for Fermi GBM and RXTE 
data respectively. 

Joint fitting of all three datasets results vcix^ed ~ ^ 

degrees of freedom and the best-fit results are presented in Fig.[2 
and Table. [T We wish to emphasize that contrary to findings by 
Zhang et al. ( |2005[ l it was possible to describe the spin evolu¬ 
tion of the source during all three outbursts without assuming 
any change in orbital parameters. Obtained results are generally 
consistent with the values reported by Zhang et al. ( 2005) 1 and 
Raichur & Paul (20101 with the exception of the orbital period 


anymore. We conclude, therefore, that there is no evidence for 
apsidal motion in the system. 


value. We note that in our case it is much better defined as the 
periastron passage time is very well constrained in both RXTE 
and Fermi datasets and the number of orbital cycles between 
the two outbursts is known. Other discrepancies likely arise due 
to the different assumptions on the intrinsic spin frequency of 
the pulsar coupled with insufficient orbital phase coverage and 
limited accuracy of individual spin frequency measurements in 
previous studies. 


Implications for the binary orbit. We have shown that it is pos¬ 
sible to describe the data outbursts in 1983, 2005 and 2015 as¬ 
suming no change in the orbital parameters, thus resolving the 
discrepancy in a sin(/) values reported by Stella et akll 1985) 1 and 
Zhang et al. l |2005 1. The discussion by [Zhang et al. ( 2005| l regard- 
ing a possible apsidal motion in the system is thus not required 


On the other hand, the low a sin i value reported by Stella 
et al. (1985 i led Negueruela et al. ( 1999|l to conclude that the 


Be companion in V 0332+53 must be strongly undermassive. 
Indeed, assuming asint ~ 50Its the projected rotational ve¬ 
locity of the Be star measured by [Negueruela et al. (1999 1 
Vt sin/t ~ 100 - 200 km s * implies t, > 12° - 24° if the in¬ 
trinsic rotational velocity of the star is below the break-up value 
V, < 0.8t>break- up ~ 4 80kms~*. For orbital parameters reported 
by [Stella et al.[ ( 1985) l this indeed implies a strongly undermassive 
optical companion with <5-7 Mq unless its equatorial plane 
is misaligned with the orbital plane of the system ( [Negueruela 
et al.|[1999jl. For larger a sin / ~ 78 It s found in this work the 


same considerations considerations imply < 8 - 50 Mq, i.e. 
compatible with the expected mass of a Be star. We conclude, 
therefore, that there is also no evidence for a tilt between its 
equatorial plane and orbital plane of the system. 

Intrinsic spin evolution Once the orbital parameters of the sys¬ 
tem are determined it is possible to assess the intrinsic spin evolu¬ 
tion of the pulsar and thus probe the accretion torques acting onto 
the neutron star. To estimate the time spin frequency derivative 
we correct the observed frequency values for motion in the binary 
system using the ephemeris obtained above and calculate the fre¬ 
quency derivative comparing the measurements in adjacent time 
intervals (propagating the uncertanties). The average spin-up rate 
of Pspin ~ 5 X 10“® s d“* is comparable with Pspm ~ 8 x 10“® s d“' 
reported by Zhang et al. ( 2005| l and, as discussed by these authors, 
is in rough agreement with the accretion torque theory. 

The spin-up rate of the neutron star is expected to be cor¬ 
related with the accretion rate. Therefore, we estimate also the 
accretion rate for each time interval where the spin-up rate was de¬ 
termined. In particular we multiply the average Swift BAT count 
rate in 15-50keV range during respective time interval by fac¬ 
tor of 1.6 X lO^^erg s *cm“^ to calculate the flux. This factor 
was calculated based on the comparison of the Swift BAT count 
rates and flux in 3-lOOkeV range measured during the pointed 
NuSTAR observations on MJD SllXi, 57275 and 57281 (detailed 
analysis of pointed observations will be presented elsewhere). 
The luminosity and accretion rate can then be estimated assum¬ 
ing the distance to the source of 7 kpc and standard neutron star 
parameters. We find that the spin-up rate is indeed correlated with 
the observed X-ray flux as illustrated in Fig. 

Several accretion torque models describing the spin evolution 
of a neutron star accreting from a disk have been proposed. The 
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Parameter 

This work 

RPIO 

Z05 

Porb, d 

33.850(1) 

36.5(3) 

34.7(4) 

a sin(i), It s 

77.81(7) 

82.5(9) 

86(10) 

e 

0.3713(8) 

0.417(7) 

0.37(12) 

oj, deg 

277.43(5) 

283.5(9) 

283(14) 

PpA, mjd 

57157.88(3) 

53330.58(6) 

53367(1) 

AT, km 

53.97(5) 

- 

59(7) 

f(M), Mo 

0.441(1) 

- 

0.58(23) 


Table 1. The best-fit orbital parameters of V 0332+53 including 
the Fermi GBM data from 2015 outburst (left column). Values 


obtained by Raichur & Paul (20101 (RPIO) and Zhang et al. 
( |2005| ) (Z05) are also shown for reference. All uncertainties are 
at Icr conbdence level. 


xl0“"^ 



Fig. 2. Observed spin-up rate as function of source X-ray lumi¬ 
nosity for the rising (circles) and declining (diamonds) parts of 
the 2015 outburst. Model predictions after Ghosh & Lamb ( 1979^ 
and|Parfrey et al.|(|2015|) are also shown for reference. 


model by Ghosh & Lamb (19791 invokes angular momentum 
transport from a neutron star to an accretion disk threaded by the 
magnetic field lines. As can be seen in Fig.|^ this model predicts 
somewhat higher spin-up rate than observed (assuming the mag¬ 
netic field of the neutron star of B - 3.5 x 10*^ G corresponding 
to the highest cyclotron line energies reported by|Tsygankov et alT 


2006|l. The discrepancy is not large, however, and could be rec¬ 


onciled if we assume that the accretion rate is overestimated by 
factor of two. 


In a more recent study |Parfrey et al.1 ( |20l5| l argue that the an¬ 
gular momentum is removed from the neutron star by the electro¬ 
magnetic outflow along the field lines opened by differential rota¬ 
tion of the magnetic field and the accretion disk. The total torque 
in this model depends strongly on the magnetospheric radius, 
which is assumed to constitute a fraction of the Alfvenic radius 
fa = (/r"^/2GMM^)F’ for spherical accretion, i.e. The 

value of ^ ~ 0.5 is expected from MHD simulations (Long et al. 
|2005||Bessolaz et al.|2008||Zanni & Ferreira|2013| l. However, in 
this case the model predicts significantly higher spin-up rate for 
V 0332+53 than observed. To match the observations with the 
model prediction one can assume either that the accretion rate is 
over-estimated by factor of ten (which is unlikely) or that the ^ 
must be reduced by factor of two (i.e. f = 0.25 as presented in 
Fig.|^ which could signify that the disk pushes further into the 
magnetosphere than expected. 


3. Conclusions 


In this note we analyzed the Fermi GBM spin history of 
V 0332+53 during the giant outburst in 2015 together with histor¬ 
ical data from previous outbursts in 1983-1984 and 2004-2005. 
Our results are generally consistent with earlier estimates by 
Zhang et al. ( 2005| l; Raichur & Paul 


(20101 and similar to the 


preliminary findings by the Fermi GBM team (based only on the 
data from the latest outburst). 

For the first time we succeeded to describe the spin evolu¬ 
tion of the source during all three outbursts with no change in 
orbital parameters between them thus resolving a long standing 
discrepancy between the orbital solutions reported by Stella et al. 
( |1985 I, [Zhang et ah ( 2005| l, and Raichur & Paul (|2010| l. Our 
finding makes the discussion by Zhang et al. ( |20()5 I regarding 
the possible apsidal motion in the system unnecessary. We note 
also that in the light of updated ephemeris the conclusion by 
Negueruela et al. ( 1999| l regarding the misalignment of orbital 
plane of the system and Be stars equatorial plane does not hold 
anymore. We were also able to significanlty improve the accuracy 
of the orbital solution. 

We find the the intrinsic spin evolution of the pulsar to be 
complicated with the spin-up rate being correlated with the ac¬ 
cretion rate. The observed spin-up is qualitatively consistent with 
existing torque models assuming the neutron star has the mag¬ 
netic field of ~ 3 X 10^^ G as determined from cyclotron line 
energy although model uncertanties and the uncertainty in accre¬ 
tion rate of the source prevent us from any conclusions regarding 
preferred torque model. 
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